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Abstract 

The density matrix renormahzation group method is used to investigate the Peierls transition 
for the extended Hubbard model coupled to quantized phonons. Following our earlier work on 
spin-Peierls systems, we use a phonon spectrum that interpolates between a gapped, dispersionless 
(Einstein) limit to a gapless, dispersive (Debye) limit to investigate the entire frequency range. 
A variety of theoretical probes are used to determine the quantum phase transition, including 
energy gap crossing, a finite size scaling analysis, and bipartite quantum entanglement. All these 
probes indicate that a transition of Berezinskii-Kosterlitz-Thouless-type is observed at a non-zero 
electron-phonon coupling, g^, for a non- vanishing electron-electron interaction. 

An extrapolation from the Einstein limit to the Debye limit is accompanied by an increase in 
gc for a fixed optical {q = vr) phonon gap. We therefore conclude that the dimerized ground 
state is more unstable with respect to Debye phonons, with the introduction of phonon dispersion 
renormalizing the effective electron-lattice coupling for the Peierls-active mode. 

By varying the Coulomb interaction, U , we observe a generalized Peierls transition, intermediate 
to the uncorrelated {U = 0) and spin-Peierls {U — t- co) limits, where U is the Hubbard Coulomb 
parameter. 

Using the extended Hubbard model with Debye phonons, we investigate the Peierls transition 
in irans-polyacetylene and show that the transition is close to the critical regime. 

PACS numbers: 71.10.Fd, 71.30. +h, 71.38.-k 
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I. INTRODUCTION 



Low- dimensional electronic materials are known to be highly susceptible to electron- 
phonon-driven structural distortions. Almost half a century ago, Peierls demonstrated that 
a one-dimensional (ID) metallic system can support a periodic modulation in the equilibrium 
positions of the lattice ions-. For the case of a half-filled band, the broken-symmetry phase 
is commensurate with the lattice, resulting in a doubling of the unit cell of the ground 
state (GS). Since the dimerization opens a gap at the Fermi surface, the Peierls process 
transforms the metal to a dielectric phase, with the increase in lattice energy associated 
with the permanent distortion being offset by the reduction in electronic kinetic energy. 
Spontaneous dimerization has been noted in many quasi- ID materials, ranging from organic 
conjugated polymers-!^ and charge-transfer salts^ to inorganic blue bronzes^ and MX-chains^. 

The Peierls instability is well understood in the static lattice limit for which the frequency, 
Uj^, of the Peierls-active mode is taken to be much smaller than the electron hopping in- 
tegral, t. In the adiabatic phonon limit, the GS is known to have a broken-symmetry 
staggered dimerization for arbitrary electron-phonon coupling. Experimentally, such be- 
havior was first observed in the 1970s for the organic compounds of the TTF and TCNQ 
series^. For many quasi-one-dimensional materials, however, the zero-point fluctuations of 
the phonon field are comparable to the amplitude of the Peierls distortion^"^. Lattice dy- 
namic (quantum phonons) effects should therefore be included in a full theoretical treatment 
of Peierls-distorted systems. 

Likewise, interest in models of spins coupled dynamically to phonons increased signifi- 
cantly when it was shown that the first inorganic spin-Peierls (SP) compound CuGeO^ 
exhibits no clear scale-separation between magnetic and phononic energies. Moreover, in 
contrast to the organic SP materials, no phonon-softening is observed at the transition. SP 
physics, then, is plainly in the non-adiabatic regime. 

Using the density-matrix renormalization group (DMRG), it has been demonstrated 
that quantum fluctuations destroy the Peierls state for small, non-zero couplings in both 
the spinlessii and spin-i^ Holstein models at half- filling. Analogous results for the XY- 
SP model with gapped, dispersionless (Einstein) phonons were obtained by Caron and 
Moukouri^-, using finite-size scaling analysis of the spin gap to demonstrate a power-law 
relating the critical coupling and the Peierls-active phonon frequency: g^"^ ~ wJJ'^. Citro et 
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alM used a renormalization group (RG) treatment of the bosonized Heisenberg-SP model to 
demonstrate similar behavior in the antiadiabatic phonon regime {utt/J » 1). In general, 
for models with sufficiently large Einstein frequency, gapped phonon degrees of freedom can 
always be integrated away to generate a low-energy effective-fermion Hamiltonian charac- 
terized by instantaneous, non-local interactions^ii^. 

Gapless, dispersive (Debye) phonons were found by the present authors to destabilize the 
broken-symmetry GS of the Heisenberg-SP chain^^, which maps to a spinless-fermion-Peierls 
chain under Jordan- Wigner (JW) transformation. By interpolating between the Einstein- 
and Debye-phonon limits, the spin-Peierls phase was shown to be more unstable with respect 
to dispersive lattice degrees of freedom. For the Su-Schrieffer-Heeger (SSH) model, Fradkin 
and Hirsch undertook an extensive study of spin-i (n = 2) and spinless (n = 1) fermions 
using world-line Monte Carlo simulations^^. In the antiadiabatic limit (i.e. vanishing ionic 
mass M), they mapped the system onto an n-component Gross-Neveu model, known to 
exhibit long-ranged dimerization for arbitrary coupling for n > 2 (although not for n = 1). 
For M > an RG analysis indicates the low-energy behavior of the n = 2 model to be 
governed by the zero-mass limit of the theory, indicating that the spinful model presents a 
dimerized GS for arbitrarily weak e-ph couplings^^. 

Although of theoretical interest, independent-electron models are not sufficient to give a 
quantitative account of the properties of physical systems exhibiting a Peierls distortion - 
for that, electron-electron (e-e) interactions must to included. The interplay between e-e and 
electron-phonon (e-ph) interactions results in an extremely rich phase-diagram of broken- 
symmetry GSs, each supporting a range of low-energy electron-lattice excitations: solitons, 
polarons, lattice "breathers", eto^>^. For the half-filled Hubbard model, repulsive on-site 
interactions in the charge sector {U > 0) give rise to the opening of a charge gap A*^'^^ and, 
in the absence of e-ph couplings, the system is a Mott insulator (MI): a critical dielectric 
phase exhibiting algebraically decaying spin-spin correlations^^. Hence, the transition from 
the Mott insulator (MI) to the Peierls insulator (PI) phase is accompanied only by the 
generation of a spin gap A*^'^-', the charge gap having arisen by virtue of mutual electronic 
repulsions. Longer range e-e interactions, e.g. next-nearest neighbor terms {V), can destroy 
the Mott state, however, and are expected to influence the MI-PI transition. Even in the 
absence of lattice degrees of freedom, the phase diagram of the simplest half-flUed extended 
Hubbard model is still controversial: its behavior close to the U = 2V line has attracted 
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significant attention, with a Peierls-like bond-ordered GS predicted to exist^. 

Sengupta et al.—, using the extended Hubbard model coupled to Einstein bonds phonons, 
demonstrated the destruction of the PI below a non-zero value of the e-ph coupling, in 
agreement with earlier independent-electron treatments. Work by Zimanyi et al?^ on one- 
dimensional models with both e-e and e-ph interactions indicated the development of a 
spin gap provided the combined backscattering amplitude gj = giiuj) + Qiioj) < 0, where 
gi{(jj) is the contribution from e-e interactions and ^i(w) < is the e-ph contribution in the 
notation o^. Hence, for the pure spinful SSH model {U = V = 0) , gi = and gj<0 for 
any nonzero e-ph coupling, implying a Peierls GS for arbitrary e-ph coupling, in agreement 
with the earlier MC results^. It should be noted, however, that a dimerized state was also 
predicted for g^ = O"*" in the spin-i Holstein model, for which later large-scale calculations 
indicated a nonzero critical coupling2^. 

In this paper we examine the influence of gapless, dispersive phonons on the GS of the 
extended Hubbard-Peierls (EHP) chain. We explicitly probe the MI-PI transition, studying 
the model as a function of the electron-phonon interaction g, Coulomb interactions U and 
V, and phonon frequency Utt- In all cases we consider a half-filled band. Since the parameter 
space of this model is rather large, we restrict our study to a physically reasonable ratio 
U/V = 4 of the e-e parameters for the full phonon-frequency range. In addition, we investi- 
gate the model for w^/ 1 = 1 for all values of U. The f/ — t- oo limit is of particular relevance 
to our earlier work since^^, in the limit of large Hubbard interactions, the EHP Hamiltonian 
maps onto the quantum Heisenberg- Peierls antiferromagnet. However, for most polyenes, 
the Coulomb interactions are not large enough to justify the spin model^^. We also note that 
in the limit of vanishing phonon frequency (M — )■ oo) the model maps onto the extensively 
studied classical adiabatic chain. 

That the dispersive-phonon EHP model is yet to receive the same level of attention as its 
gapped, dispersionless counterpart is due in part to the presence of hydrodynamic modes, 
resulting in logarithmically increasing vibrational amplitudes with chain length. To this 
end, acoustic phonons have been assumed to decouple from the low-energy electronic states 
involved in the Peierls instability, motivating the retention of only the optical phonons close 
to g = TT— In this regard, optical phonons have been expected to be equivalent to 
fully quantum mechanical SSH phonons. Even for pure Einstein phonons, however, Wellein, 
Fehske, and Kamp^ found that the singlet-triplet excitation to be strongly renormalized 
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when phonons of all wavenumber are taken into account, the restriction to solely the q = tt 
modes leading to a substantial overestimation of the spin gap. Physically, this implies that 
the spin-triplet excitation is accompanied by a local distortion of the lattice, necessitating 
a multiphonon mode treatment of the ionic degrees of freedom. Our recent work on the SP 
model has indicated that truncating the Debye-phonon spectrum, leaving only those modes 
which couple directly to the Peierls phase, is not physically quantitatively reasonable. 

We use the DMRG^ technique to numerically solve the EHP model for t = U/A = V 
with a generalized gapped, dispersive phonon spectrum. The phonon spectrum interpolates 
between a gapped, dispersionless (Einstein) limit and a gapless, dispersive (Debye) limit. 
We proceed by considering a system of tightly-bound Wannier electrons dressed with pure 
Einstein phonons for which we observe a Berezinskii-Kosterlitz-Thouless (BKT) quantum 
phase transition at a non-zero electron-lattice coupling. Progressively increasing the Debye 
character of the phonon dispersion (at given optical phonon adiabaticity) results in an in- 
crease in the critical value of the e-ph coupling, with the transition remaining in the BKT 
universality class (see Section [III CI) . These findings are corroborated by an array of inde- 
pendent verifications: energy-gap crossings in the spin excitation spectra (see Section UlIAl) . 
finite-size scaling of the spin-gap (see Section FlIIBI) . and quantum bipartite entanglement 
(see Section UlI pp . Our approach here is equivalent to that described in ref^^ to investigate 
the BKT transition in the spin-Peierls model. 

We note that earlier DMRG investigations of the EHP Hamiltonian with Debye phonons 
indicated a dimerized GS for arbitrary coupling^- . This conclusion was based on the behavior 
of the staggered phonon order parameter, which we have since shown to be an unreliable 
signature of the transitionii. We therefore pursue alternative characterizations of the Peierls 
state in this work. 

In the next Section we describe the model, before discussing our results in Section lllli 

II. THE MODEL 

The extended Hubbard-Peierls Hamiltonian is defined by, 

H = He-e + He-ph + -f^ph- (l) 
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ife-e describes the electronic degrees of freedom, 



l,cr 

+UJ2 (ni, - (^Nij, - +Vj2iNi - l){Ni^i - 1), 



and -ffe-ph the e-ph couphng, 



(2) 



l,a 



Here, Ni^j = cj^cia, where c| (qo-) creates (annihilates) a spin-cr electron at Wannier site / of 
an A^-site ID lattice, ui is the displacement of the Zth ion from equilibrium, a is the e-ph 
coupling parameter, and U and V are the on-site and nearest neighbor Coulomb interactions, 
respectively. 

ifph describes the lattice degrees of freedom. In the Einstein model the ions are decoupled, 

(3) 
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In the Debye model, however, the ions are coupled to nearest neighbors. 
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For the Einstein phonons it is convenient to introduce phonon creation, 6j, and annihila- 
tion operators, 6/, for the Zth site via, 

1/2 



h 



2Mu 



X 



(bl + k 



and 



where 



Making these substitutions in Eq. (|2]) and Eq. ([3]) gives. 
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(c^Q+ia + H.c.) 



(8) 
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and 



(9) 



where Bi = |(&|^ + 6;) is the dimensionless phonon displacement and 
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(10) 



is the dimensionless e-ph coupling parameter. 

For the Debye phonons we introduce phonon creation and annihilation operators defined 
by Eq. ([5]) and Eq. ([6]) where 



Making these substitutions in Eq. (|2]) and Eq. (jl]) gives, 



H, 



c-ph 



(cLc/+i<7 + H.c.) 



and 



where, 



H^^ = nwoJ2 i^i^^ + 1)- E 
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may be diagonalized by a Bogoluibov transformation^^ to yield. 



(11) 



(12) 



(13) 



(14) 



(15) 



where ujo^q) is the dispersive, gapless phonon spectrum. 



UD{q) = V2ujd sin ^ 



(16) 



for phonons of wavevector q. 

We now introduce a generalized electron-phonon model with a dispersive, gapped phonon 
spectrum, via 



He-ph — t E^ 



1 + 9 



t 



{Bi — Bi^i 



(17) 
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and 

i^ph = Hue + cod) J2 (^f^' + ^U^h (18) 

Again, Eq. ( ITS]) may be diagonalized to give, 

ifp = /i^^u;(g)/3^/3g + constant, (19) 

where, 

u{q) = (w£; + Wd) ( 1 - ( — — — ) COS g ) , (20) 

is the generahzed phonon dispersion, as shown in Fig. [H 
The g = phonon gap frequency is, 

uj{q = 0) = Wo = ij^Eij^E + ^d))^^"^ (21) 
and the q = tt optical phonon frequency is, 

1 /2 

uj{q = tt) = ui^ = {{ue + ujD)i'2uE + ujd)) ■ (22) 
We now define the dispersion parameter 7 as, 

7 = Wo/w^- (23) 

7 is a mathematical device that interpolates the generalized model between the Einstein 
(7 = 1) and Debye (7 = 0) limits for a fixed value of the q = it phonon frequency, Ut^. The 
dimensionless spin-phonon coupling, g, as well as ooT^/t and 7 are the independent parameters 
in this model, ue and ud, on the other hand, are determined by Eq. (!2T1) . ( l22l) . and (!23l) . 

The generalized model can be mapped onto the physical Einstein and Debye models by 
the observation that in the Einstein limit, 

u]^ = ue = ujo] (24) 
9 = 9e, 

while in the Debye limit, 

cJtt = \/2ujd = 2ci;o; (25) 
9 = 9d/2'/\ 
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q/n 



FIG. 1: Generalized phonon dispersion, defined in Eq. ()20p . (1 — 'y)uT^ is the phonon 'band width' 
(which vanishes in the Einstein-hmit), while "jlVt^ is the phonon 'mass-gap' (which vanishes in the 
Debye-limit). 

The introduction of a generalized phonon Hamiltonian avoids the problems associated 
with hydrodynamic modes and places a criterion on the reliability of the gap-crossing char- 
acterization of the critical coupling (as described in Section [III A|) . Starting from the ex- 
tended Hubbard-Peierls Hamiltonian in the Einstein limit (7 = 1), the effect of dispersive 
lattice fluctuations can be investigated via a variation of 7. The Debye limit is then found 
via an extrapolation of 7 — ?■ 0. 

We note that the EHP model is invariant under the particle-hole transformation, c]^ — )■ 
(— l)Vj5-. This so-called charge-conjugation symmetry is exact for vr-electron models but it is 
only an approximate symmetry for conjugated polymers and is strongly violated for systems 
possessing heteroatoms. Nevertheless, for the EHP model at half-filling it is expedient to 
employ particle-hole symmetry to distinguish between different types of singlet excitations, 
as described in §111 A[ 

The many-body problem is solved using the density matrix renormalization group 
(DMRG) method^ with periodic boundary conditions throughout. Our implementation 
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of the DMRG method, including a description of the adaptation of the electron-phonon 
basis and convergence, is outhned in^iili^S. 

III. RESULTS AND DISCUSSION 
A. Gap-crossing 

Since both the MI and PI possess a non- vanishing charge gap, A'^'^^ spectroscopy of the 
spin excitation sector must be used to characterize the GS phase. For the Einstein model 
with a non- vanishing value ofuE, the critical e-ph coupling, Qc, may be determined using the 
gap-crossing method of Okamoto and Nomura— (as illustrated in Fig. 2 of ref^^). If the iV-site 
system is a MI exhibiting quasi- long-range order for < g < gdN), the lowest spin-sector 
excitation is to a triplet state, i.e. A^t < Ass and limAr^oo ^st = limTv^oo Ass = 0, where Ast 
and Ass are the triplet and singlet gaps, respectively. Conversely, for g > gd^N), the system 
is dimerized with a doubly-degenerate singlet GS in the asymptotic limit (corresponding to 
the translationally equivalent 'A' and 'B' phases), while the lowest energy triplet excitation 
is gapped. However, for finite systems the two equivalent dimerization phases mix via 
quantum tunneling, and now Ass < ^st, with limjv_^oo Agg = and lim7v->oo Ast = A*^*^) > 0. 
The gap-crossing condition Agt = Ass therefore defines the finite-lattice crossover coupling 
gc{N). The singlet gap with which we are concerned with here is the lowest singlet (covalent) 
excitation with the same, i.e. positive, particle-hole symmetry as the GS. Conversely, the 
charge gap, A^^\ is the lowest singlet (ionic) excitation with negative particle-hole symmetry, 
i.e. A'^'^) = Ajg. For small g > gc, A^^ = A+ < A~ . Since we are concerned only with the 
spin-excitation sector, we hereafter refer to the bulk-limit spin gap, A^*^), as A. 

For the Debye model, however, the gap-crossing method fails because of the g — ?■ 
phonons that form a gapless vibronic progression with the GS. The hybrid spectrum (shown 
in Fig. [1]) allows us to extrapolate from the pure Einstein limit to the Debye limit, as the 
lowest vibronic excitation is necessarily 70;^. Provided that Ass < u}{q = 0) = 'JUt^, the 
gap crossover method unambiguously determines the nature of the GS. We can confidently 
investigate Eq. ([1]) for (0.1 < 7 < 1) with u^^/t G [1, 10], thereby determining gc{N,'~f). A 
polynomial extrapolation of 1/A^ — )■ generates the bulk-limit critical coupling g"^ for a 
given 7 (as illustrated in Fig. 2 of re^). A subsequent polynomial extrapolation determines 
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FIG. 2: Phase diagram in the g^-w^ plane for the infinite EHP chain for the Einstein-hmit [squares]; 
extrapolation to 7 = generates the Debye-limit [circles]. V = U/4 = t. The cross (x) indicates 
the parameters relevant for trans-polyacetylene (see ^III Gp . 

the 7 = (Debye) limit. A phase diagram for the EHP chain found in this way is shown in 
Fig. [21 Notice that for a fixed u-,^ the critical coupling is larger for the Debye model than 
for the Einstein model, showing that the quantum fluctuations from the g < vr phonons (as 
well as the q = tt phonon) destablize the Peierls state, in agreement with our earlier work 
on spin-phonon systems. 

Following Caron and Moukouri^^, as well as our more recent work on the Heisenberg-SP 
model, we propose a general power-law, relating the bulk-limit critical coupling to u-,^ for a 
given 7 of the EHP model, 

^-(u;Jt,7) = /3(7)(u;Jt)^W. (26) 

The infinite-chain values of /3, 1], and are given in Table [B We observe a non-zero 
critical coupling for all phonon-dispersion regimes 7, with the absolute value of increasing 
as 7 — )■ 0. We also find the power relation to be robust, extending well into the adiabatic 
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regime {uj-^lt « 1). 



7 


/3 


V 


9? 


1 (Einstein) 
(Debye) 


0.273 
0.332 


0.317 
0.469 


0.280 
0.331 



TABLE I: Gap-crossing determined bulk-limit values of /3(7) and rj^-j) for V = U/4: = t. The value 
of is shown for = t. See Eq. ([26|). 

We next consider the role of Coulomb repulsion by varying the on-site interaction U 
subject to U /V = 4. In the adiabatic limit, the amplitude of the bond alternation ini- 
tially increases with Coulomb repulsion. This is because electronic interactions suppress 
the quantum fluctuations between the degenerate bond-alternating phases, the alternation 
being maximized when the electronic kinetic and potential energies are approximately equal, 
namely when U ~ 4t^^. For larger Coulomb interactions, however, charge degrees of free- 
dom are effectively quenched and the EHP maps to the spin-i Heisenberg-Peierls chain with 
antiferromagnetic exchange^ J = At^/{U — V). In this regime, virtual exchange^ effec- 
tively lowers the barrier to resonance between the 'A' and 'B' phases, thereby reducing the 
dimerization and reconnecting with our earlier work— "i^. and is shown in Fig. [3] and Fig. [71 

The critical coupling, on the other hand, increases monotonically with U for both the 
Einstein and Debye-phonon limits, contrasting with the non-linear behavior of the spin gap. 
The f/ = intercept is found to be zero in both indicated in Fig. HI in agreement 

with earlier work on the half-filled spinful SSH model^"^. The authors of^^ tentatively 
propose power-law scaling for the critical coupling g'^ ~ U*^'^ ioi V = U/A = t and = t, 
but cite the smallness of the spin gap as (f/, y) — ?■ as a source of uncertainty below 
U = OAt. We note a similar power-law behaviour, 

^?r(f/A,7)=/3(7)(f/A)^^"\ (27) 

which implies = Vcjvr when f/ = in agreement with the earlier findings of Fradkin 
and Hirschp^ and the many-body valence bond treatment of Dixit and Mazumdar-^. The 
infinite-chain values of /3 and fj are given in Table Ull 



12 



5.5 

5 
4.5 

4 

3.5* 

3 
2.5 

2 
1.5 

10 

0.5 




o 



o 



o 



o 



o 



o 



C) 







8 

U/t 



10 12 14 



16 



FIG. 3: Bulk spin gap, A, versus U/t for the infinite EHP chain for 7 = 1 (Einstein) [squares]; 
extrapolation to 7 = generates the Debye-limit [circles]. A is evaluated for 5 = 4, i.e. well into 
the Peierls phase. lUt^ = V = U/A = t. 



7 


/3 


V 


1 (Einstein) 
(Debye) 


0.171 
0.195 


0.281 
0.401 



TABLE II: Gap-crossing determined bulk-limit values of and r/(7) for uj-^ = V = U/i = t. See 
Eq. (HI 



B. Finite-size scaling 



In order to ascertain the analytic behavior of the spin gap from the numerical data 
it is necessary to account for finite-size effects. We assume that the (singlet-triplet) gap 
An = Ast for a finite system of sites obeys the finite-size scaling hypothesis^"^ 
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FIG. 4: Phase diagram in the g^-U plane for the infinite EHP chain for 7 = 1 (Einstein) [squares]; 
extrapolation to 7 = generates the Debye-limit [circles], u.,^ = V = U/4: = t. 

with Aoo the spin-gap in the bulk limit. Recalling that = limN^oo gc{N), it follows that 
Aoo(fi'^) = and so curves of NAjy versus g are expected to coincide at the critical point 
where the bulk-limit spin-gap vanishes, as confirmed in Fig. O 

The finite-size scaling method is more robust than the gap-crossing approach, being 
applicable to the EHP Hamiltonian for all values of 7. On the other hand, its use as a 
quantitative method is limited by the accuracy with which plots may be fitted to Eq. fl28l) . 
In practice, plots of NAst{N) versus g become progressively more kinked about the critical 
point as 7 — )■ 0. Nevertheless, we find F to be well approximated by a rational function and 
the resulting (7^(7) to be in accord with the predictions of the gap-crossover method. 



C. Berezinskii-Kosterlitz-Thouless transition 



For a BKT transition the spin-gap A = limiv_^oo Ast is expected to exhibit an essential 
singularity at g^ with plots of Ast versus g for N ^ 00 found to be well fitted by the Baxter 
formS (as shown in Fig. [6]), 

A^af{g)exp{-b[f{g)]') (29) 
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FIG. 5: A^Ast(A^) versus the e-ph coupling, g, for the 7 = (Debye) EHP model for = 16 
(diamonds), 40 (squares), 100 (triangles), and 160 (circles). The curves converge at Qc (the value 
shown is obtained via gap-crossing), uj^^ = V = U/4 = t. 

whereii, 

f{9)^{9-9n-'^'- (30) 

Extrapolating Ast(A^) for 1/A^ — > generates A for a given 7 and it is possible, in 
principle, to distinguish MI from PI GSs by examining the scaling behavior of Ast(A^), 
which tends to zero in the bulk-limit for the critical MI and to a non-zero A for the gapped 
dimerized phase. However, not only must three parameters (a, b, and g"^) be obtained 
from a non-linear fit (shown in Table [In|) . but there is considerable difficulty in determining 
A accurately near the critical point: the spin-gap is extremely small even for values of g 
substantially higher than g^ due to the essential singularity in Eq. fj29|) . Determining such 
small gaps from finite-size scaling is highly problematic, with very large lattices required to 
observe the crossover from the initial algebraic scaling (in the critical spin-sector regime) to 
exponential scaling (for gapped systems). Hence, the gap-crossover method is expected to 
be substantially more accurate than a fitting procedure for the determination of the critical 
coupling, the latter tending to overestimate g^'^, as confirmed by a comparison of Tables |T] 
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FIG. 6: Bulk-limit singlet-triplet gap, A, as a function of the e-ph coupling, g, with 7 = 1 (Einstein) 
[squares] and 7 = (Debye) [circles] for = F = C//4 = t. Plots are fitted to the BKT form (Eq. 

(USD). 

and [ml 



7 


a 


b 


5c°° 


1 (Einstein) 
(Debye) 


18.501 
6.704 


2.521 
2.091 


0.285 
0.349 



TABLE III: Baxter-equation parameters obtained by fits to Eq. (I29p for o;^ = y = [7/4 = 

Finally, we consider the effect of Coulomb repulsions on the Baxter-equaton parameters 
for which the corresponding plots are shown in Fig. [71 We note the greatest amplitude. A, 
arises for U ~ At, in agreement with Section [III A[ For strong coupling, i.e. U/t = 16, the 
function approximates to that of the previously considered spin-Peierls modeU^. Beyond 
U/t ^ 16, the spinless fermion picture becomes increasingly appropriate, in agreement with 
findings for the Hubbard model^^. 
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FIG. 7: Bulk- limit singlet-triplet gap, A, as a function of e-ph coupling, g, for 7 = 1 (Einstein), 
oJtt = y = t. U/t = (triangles), U = t (circles), U = At (squares), U = IQt (diamonds), and the 
Baxter-equation fitting for the spin-Peierls casei^ {U/t ^ 00) (dashed line) are shown. 

D. Quantum bipartite entanglement 

Entanglement has been shown to play an important role in the quantum phase transitions 
(QPT) of interacting lattices. At the critical point — as in a conventional thermal phase 
transition — long-range fluctuations pervade the system. However, because the system is at 
T = 0, the non-degenerate GS is necessarily pure, indicating that the onset of (long-range) 
correlations is due to scale-invariant entanglement in the GS. 

For an A^-site lattice, bipartite entanglement is quantified through the von Neumann 
entropy^, 

Sl = -TrgpsiL) log2 psiL) = - ^i^^ hg^ i^a, (31) 

a 

where ps{L) is the reduced-density matrix of an L-site block (typically coupled to an L-site 
environment S such that 2L = N) and the z/q, are the eigenvalues of ps{L). Provided the 
entanglement is not too great and the Va decay rapidly, a matrix-product state is then a 
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good approximation to the GS^. 

Wu et al— argued, quite generally, that QPTs are signalled by a discontinuity in some 
entanglement measure of the infinite quantum system. In^ we demonstrated that the non- 
critical (gapless) phase entanglement is characterized by the saturation of the von Neumann 
entropy with increasing L, Sl growing monotonically until it saturated for some block length 
Lq, in agreement with^. The critical (gapped) phase, on the other hand, was found to ex- 
hibit logarithmic divergence in Sl at large L. 

The principal difference between qubit lattices (e.g spin chains) and itinerant-particle sys- 
tems is the requirement of wavefunction anti-symmetrization for indistinguishable fermions, 
which implies a Hilbert space lacking a direct-product structure. Such a structure may be 
recovered, however, by passing to the occupation number representation of local fermionic 
modes^: the A^-site lattice is spanned by the 4^ states {(^^ l'^);}, where site / has local 
basis {1?^)/} = {|0)/, I t)/) I i)h I ti)/}- Under this formahsm, the von Neumann entropy for 
pure states remains a well-defined entanglement measure, having been used to determine 
the phase diagram of the extended Hubbard model^. 

For a given total system size and phonon dispersion 7, the block entropy is found to be 
maximal for a non-zero spin-phonon coupling gc{N), close to the corresponding gap-crossing 
value (as shown in Table lIVp . As shown in Fig. |8l in the critical regime, g < gc{N), the 
block entropy is indeed found to scale logarithmically with system-block length, while in 
the gapped phase, g > gc{N), it is characterized by the emergence of a saturation length 
scale Lq that varies with 7. These findings are in agreement witb^ and consistent with the 
observation that the transition belongs to the BKT universality class^. 



N 


gap 

yc 


vN 

9c 


20 


0.425 


0.425 


40 


0.339 


0.339 


80 


0.335 


0.337 



TABLE IV: Consistency of the various probes of the transition: critical e-ph couplings determined 
by gap-crossing (gap) and von Neumann entropy (vN) for A^ = 20, 40, and 80 sites for the Debye 
model (7 = 0). uj^ = V = U/A = t. 
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FIG. 8: Von Neumann entropy, Sl, for the 7 = (Debye) EHP model for lattice sizes A*" = 20 
(squares), 40 (triangles), and 80 (circles); L = N/2, oj^^ = V = U/A = t. 

E. Phase diagram 

To conclude this section we discuss the phase diagram of the EHP model. Fig. |2] shows 
the phase diagram as a function of the model parameters g and the q = n phonon gap, Ut^, 
as defined in Eq. ( !T7l) and Eq. (fT9l) . Evidently, for a fixed value of the Peierls state is 
less stable to dispersive, gapless quantum lattice fluctuations than to gapped, non- dispersive 
fluctuations, implying that the q < it phonons also destablize the Peierls phase. 

It is also instructive, however, to plot the phase diagram as a function of the physical 
parameters a and Ub = K/M, as defined in Eq. (|2]), Eq. (|3]) and Eq. (|4]). The mapping 
between model and physical parameters is achieved via Eq. (fTOj) . Eq. ffT^ . Eq. (^^, and 
Eq. ( l25l) (and setting K = 1). Since u-^ = u)b for the Einstein model, whereas w,r = 2(^6 
for the Debye model, the Debye model is further into the antiadiabatic regime for a fixed 
value of Uq. We also note that for a given model electron-phonon coupling parameter, g, 
the physical electron-phonon coupling parameter, a, is larger in the Debye model than the 
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FIG. 9: Phase diagram in the a'^-Uh plane for the infinite EHP chain for 7 = 1 (Einstein) [squares] 
and 7 = (Debye) [circles]. V = U/4 = t. 

Einstein model (see Eq. f lTU]) and Eq. f lT^ ). Consequently, we expect the dimerized phase 
to be less robust to quantum fluctuations in the Debye model for fixed values of Ub and a, 
as confirmed by Fig. [91 



F. Connection to spin-Peierls systems 

Spin-Peierls chains with no net magentization, i.e. those for which Yli = 0, map to 
half-filled pseudofermion-Peierls models under a Jordon-Wigner ( JW) transformation. Using 
the renormalization group, the umklapp contribution to inter-fermion scattering is known to 
play a key role in the existence of a broken-symmetry GS, with the pseudoelectron-phonon 
coupling generating retarded backscattering (gi) and umklapp {g^ = —gi) couplings; by 
virtue of the Pauli Principle, however, the local character of the pseudofermion backscatter- 
ing cancels out^^. RG equations indicate that unless the nonlocal contribution to g^ has both 
the right sign and bare initial value, the umklapp processes are irrelevant and the quantum 
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system gapless. Conversely, if the threshold condition is satisfied, the umklapp processes 
and vertex function grow to infinity, signaling the onset of gapped excitations. 

The spin gap in spinful fermionic systems, on the other hand, arises because of attractive 
overall backscattering, with the sign of gj determining the existence (nonexistence) of an 
electronic Peierls GS--. The result of the commensurability effects arising from the relevance 
of the umklapp term for the half-filled band is the concurrent opening of a charge gap A^'^\ 
separating the well-known Hubbard sub-bands. Coupling to quantized phonon degrees of 
freedom renormalizes both backward and umklapp terms. 

Recapitulating the original treatment of Peierls for the half-filled bandi, U < 2V favors 
singly-occupied lattice sites. In the Luttinger liquid phase, then, we have one electron per 
unit cell and lattice constant a = TT/2kp. For "small" U (and taking g to be critical), 
coupling to a distortion of wavevector 2fci? causes a spin gap to open spontaneously. The 
unit cell doubles in size a — )■ a' = 2a, accommodating two electrons. We have, then, 
that a' = 27T /2kF = ir/kp, i.e. the reciprocal lattice vector and Fermi wavevector are 
coincident, opening a gap at the Fermi surface. Proceeding to the limit U = oo, double 
occupancy is prohibited, which decouples spin and charge dynamics, effectively quenching the 
charge degrees of freedom and giving rise to an essentially filled valence band: k^^°° = n/a 
and hence kp^°° = 2kF- The system, under a JW transformation, maps to a spinless 
tight-binding fermion chain, i.e. one pseudoelectron per unit cell, the mapping generating 
an alternating real-space occupancy pattern reminscent of the Akp charge-ordered state in 
quarter-filled spin-1/2 systems. However, since the undistorted chain has one pseudo fermion 
per unit cell, the system is half-filled and umklapp scattering becomes relevant above a 
certain pseudoelectron-phonon-coupling threshold, opening a mass gap in the spectrum. 
Mapping the system back to a spin chain under an inverse JW transformation corresponds 
to a spin- gapped phase. In this way, we have a unified treatment of the electronic Peierls 
and spin-Peierls phases. 

G. frans-polyacetylene 

Electron-lattice and inter-electron interactions in tt conjugated systems, such as trans- 
polyacetylene (t-PA), are conveniently modelled by the EHP chain for 7 = 0, i.e., the 
extended Hubbard-SSH system (EH-SSH)^>^'^. The low-energy electronic properties are 
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dominated by a single, half-filled band involving the C2p^ orbitals. vr electrons, interacting 
via long-range Coulomb forces, are coupled to longitudinal phonons, with changes in bond 
length leading to linear corrections to the hybridization integrals. The interplay between the 
delocalization of the valence electrons and the associated local fluctuations of the Coulomb 
repulsion energy is fundamental in determining the dimerization of t-PA, which has been 
successfully described as a Mott-Peierls system by the EHP model^*^. 

Dimerization in t-PA has also been studied in the adiabatic limit with the Pariser-Parr- 
Pople-Peierls model using the model parameters t = 2.539 eV, U = 10.06 eV, uq = 0.2 eV for 
C-C stretches, and the electron-phonon parameter A = 2a'^/7[Kt = 0.115^'^^. The relevant 
parameters for the EH-SSH model are then u.,^ = 0.158t, V ^ U/4 ^ t, and (using Eq. (IT^ 
and Eq. fl2^ ) gt-pA = 0.253. The critical value of g for these parameters is gc = 0.140. 

The cross on Fig. |2] indicates t-PA's position in the phase diagram using these parameters. 
Fig. [TOl shows the spin gap versus an arbitrary value of g (with other parameters fixed) and 
an arbitrary value of (with other parameters fixed). Evidently, although the value of the 
bulk spin gap is close to its asymptotic value as function of g, t-PA is close to the critical 
regime, resulting in large quantum fluctuations of the bond lengths^i"-. 

IV. CONCLUSIONS 

The interplay between the repulsive, instantaneous Coulomb interactions and the at- 
tractive, retarded interactions mediated by phonons in a ID tight-binding electron system 
results in competition between the Mott insulator and Peierls insulator ground states. For 
the extended Hubbard-Peierls chain the former becomes unstable with respect to lattice 
dimerization above a non-zero e-ph coupling threshold for all phonon gaps, 'ju^^. This ob- 
servation is true for antiadiabatic phonons (t/^Tr << 1) and remains applicable well into the 
adiabatic region of phonon phase space {uT^/t < 1). 

Increasing the contribution of dispersive phonons to ifph for fixed Coulomb interaction 
gives rise to a monotonic increase in the critical coupling for all Coulomb repulsions U, sup- 
porting the intuition that gapless phonons more readily penetrate the ground state (with the 
q < TT phonon modes renormalizing the dispersion at the Peierls- active modes). This obser- 
vation has been corroborated by an array of independent verifications and is in agreement 
with our previous work on the spin Peierls chain^. 
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FIG. 10: Bulk-limit spin gap, A, for the EH-SSH model of t-PA as a function of (a) e-ph coupling, 
g (with V = U/A = t and oj-j^ = 0.158t) and (b) phonon frequency, cj-^/t, (with V = U/4 = t and 
g = 0.253). 

The DMRG method has also been used to analyze the effect of varying U/t from the 
non-interacting hmit {U/t = 0) to the strongly correlated Heisenberg limit {U/t cxo), 
subject to U = AV. For U < 4t we observe an enhancement of the spin gap in the presence 
of repulsive interactions, with the dimerization being maximal for U approximately equal to 
the bandwidth. For larger U/t, the atomic charge fluctuations are severely reduced and the 
low-energy properties of the (quasi-localized) electrons are dominated by their spin degrees 
of freedom. The resulting spin Peierls state is regarded as arising from the alternation of 
the strength of antiferromagnetic correlations between adjacent spins. 

Finally, using the extended Hubbard model with Debye phonons, we investigated the 
Peierls transition in trans-polyacetylene and showed that the transition is close to the critical 
regime. 
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Appendix A: DMRG convergence 



We solve Eq. ( !T7|) and Eq. ( !T8|) using the real-space density matrix renormalization group 
(DMRG) metliod22, with ten oscillator levels per site, typically ~ 200 block states and ca. 10^ 
superblock states. Finite lattice sweeps are performed at target chain lengths under PBC 
In situ phonon optimization was carried out using the DMRG approach outlined in^i^. The 
convergence indicators are shown in Tables |V] and IVIl with additional convergence tables in 
ref^ for the same model. 



e 


Eg/t 


M 


SBHSS 


10- 


-10 


-54.2878527 


872 


57320 


10" 


-11 


-54.3284231 


1034 


91802 


10" 


-12 


-54.3364231 


872 


164382 


10" 


-14 


-54.3376622 


1118 


412344 


10" 


-15 


-54.3376701 


1056 


582120 



TABLE V: GS energy, Eg/t, of EHP model as a function of the density-matrix eigenvalue product 
cutoff, e, number of system block states, M, and the superblock Hilbert space size, (SBHSS) for a 
32-site chain with 10 oscillator levels per site, 7 = 1, and uJt^ = V = U/A = t. 



Number of oscillator levels per site 


N 


2 


5 


8 


10 


8 


-13.722632 


-13.816262 


-13.824165 


-13.824345 


16 


-23.381111 


-23.602836 


-23.628901 


-23.629002 


32 


-33.895366 


-34.405452 


-34.447198 


-34.447644 



TABLE VI: GS energy, Eg/t, of EHP model as a function of the number of oscillator levels per 
site for given number of sites N. e = 10~^^, 7 = 1, and lOt^ = V = U/A = t. 
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